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Abstract 

The existence of primordial adiabatic Gaussian random-phase density 
fluctuations is a generic prediction of inflation. The properties of these 
fluctuations are completely specified by their power spectrum A 2 s {k). The 
basic cosmological parameters and the primordial power spectrum together 
completely specify predictions for the cosmic microwave background radiation 
Q\ . anisotropy and large scale structure. Here we show how we can strongly 

constrain both A 2 s {k) and the cosmological parameters by combining the data 
Q-f from the Microwave Anisotropy Probe (MAP) and the galaxy redshift survey 

from the Sloan Digital Sky Survey (SDSS). We allow A 2 s {k) to be a free function, 
and thus probe features in the primordial power spectrum on all scales. If we 
assume that the cosmological parameters are known a priori and that galaxy 
bias is linear and scale- independent, and neglect non-linear redshift distortions, 
the primordial power spectrum in 20 steps in log A; to k < 0.5/iMpc -1 can be 



determined to ~ 16% accuracy for k ~ 0.01/iMpc -1 , and to ~ 1% accuracy for 
k ~ 0.1/iMpc -1 . The uncertainty in the primordial power spectrum increases by 
a factor up to 3 on small scales if we solve simultaneously for the dimensionless 
Hubble constant h, the cosmological constant A, the baryon fraction Qb, the 
reionization optical depth r ri , and the effective bias between the matter density 
field and the redshift space galaxy density field b e g. Alternately, if we restrict 
A 2 s (k) to be a power law, we find that inclusion of the SDSS data breaks 
the degeneracy between the amplitude of the power spectrum and the optical 
depth inherent in the MAP data, significantly reduces the uncertainties in the 
determination of the matter density and the cosmological constant, and allows 
a determination of the galaxy bias parameter. Thus, combining the MAP and 
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SDSS data allows the independent measurement of important cosmological 
parameters, and a measurement of the primordial power spectrum independent 
of inflationary models, giving us valuable information on physics in the early 
Universe, and providing clues to the correct inflationary model. 
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1. Introduction 

Standard cosmology poses three basic conundrums: why the universe is homogeneous 
on large scales, why the cosmological density parameter Q m is close to unity, and what 
seeds the large-scale structure in the universe. Inflation provides compelling solutions to all 
these problems ( [Kolb fc Turner 1990| ). 



The cosmic microwave background radiation (CMB) is observed to be isotropic to an 
accuracy of better than 10~ 4 . This implies that the Universe was homogeneous to this 
accuracy over many horizon scales at last scattering. Within the standard cosmology there 
is no causal process which could have created the observed homogeneity. The basic idea 
of inflation is that the observable Universe grew from an initial patch small enough to fit 
inside the horizon. It is hypothesized that there were particle physics mechanisms which 
led to the rapid expansion of the early Universe required to produce the huge growth factor 
necessary to solve the cosmological problems above. In this paper, we assume that inflation 
produces primordial adiabatic Gaussian random-phase density fluctuations due to quantum 
fluctuations in the inflaton field. The properties of these fluctuations are completely 
specified by their power spectrum ^4|(fc)0- This paper will emphasize how well the shape of 
the primordial power spectrum can be measured from CMB anisotropy and redshift survey 
data. 



There currently exists a broad range of inflationary models ( |Kolb 1997 , Turner 1997] ). 



Some of these models predict power spectra that are almost exactly scale-invariant ( |Linde| 



1983|) , or are described by a power law with spectral index less than one ([Freese, Frieman 



fc Olinto 199(J| , [La fc Steinhardt 1991] ), while others predict power spectra with slowly 



varying spectral indices ([Wang 1994]) , or with broken scale invariance ( [Holman et al.| 



1991ab , Adams, Ross, fc Sarkar 1997 , Lesgourgues, Polarski, fc Starobinsky 1997 ). 



Indeed, there are physical reasons to believe that the primordial power spectrum has 
breaks in its power- law form. In effective (i.e., not complete) theories with two scalar 
fields, inflation may occur in two stages. The two stages of inflation can generate density 
perturbations with different power-law indices, with a step in amplitude of the primordial 
power spectrum on the scale of the transition between the two phases ( [Holman et al.| 



1991ab|) . Fry & Wang (1992) found that such models can have a significant signature on 



small-scale CMB anisotropies. 

Recently, Adams, Ross, & Sarkar (1997) proposed a new multiple inflation model. Since 



3 Certain two-field inflationary models predict isocurvature as well as adiabatic 
fluctuations (|Kofman fc Linde 1987] ). Since the CMB and LSS predictions of isocurvature 
models differ significantly from adiabatic models and since these differences are not 
degenerate with parameter variation, we could, in principle, also fit for an isocurvature power 
spectrum. However, in this paper, we limit ourselves to considering adiabatic fluctuations. 
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attempts at an unified description of the strong, weak and electromagnetic interactions 
usually involve several stages of spontaneous symmetry breaking, they considered the 
effects of such symmetry breaking during an era of inflation in supergravity models. They 
found that there can be a succession of short bursts of inflation; the density perturbations 
produced during each burst is nearly scale-invariant but with differing amplitudes, and 
between each burst there is a brief period during which scale-invariance is badly broken. 

Given the range of possibilities for the primordial power spectrum, we would like to 
quantify how well it can be measured generically. Thus, we take the primordial power 
spectrum to be a free function in this paper. 

The upcoming data from the Microwave Anisotropy Probe (MAP; cf., Bennett et al. 
1997; [http : //map . gsf c . nasa . gov| ) and the Sloan Digital Sky Survey (SDSS; cf., Gunn 



& Weinberg 1996) provide a unique opportunity for constraining inflationary models. It 
is useful to expand the temperature fluctuations in the CMB into spherical harmonics: 
8T/T{i) = J2im a T,imYim(Y), where r is the unit direction vector in the sky. MAP measures 
the angular power spectrum flSeljak fc Zaldarriaga 1996] ) 



C T i = (\a TM \ 2 ) = (4vr) 2 / y A 2 s (k) \A Tl (k,r = r 



I 2 - 



where A 2 s {k) is the power spectrum of the primordial density fluctuations (defined such 
that A 2 s {k) = 1 for the scale-invariant Harrison-Zel'dovich spectrum), An(fc,r = r ) is an 
integral over conformal time r of the sources which generate the CMB fluctuations, and To 
is the conformal time today. 

The power spectrum of mass fluctuations in the linear regime today is 

P(k)=P kA 2 s (k)T 2 (k), (2) 

where Pq is a normalization constant, and T(k) is the transfer function, which depends on 
physics at matter-radiation equality and decoupling. Note that Pq has been absorbed into 
the definition of ATi(k, r = r ) in Eq. ([[]). 

The galaxy redshift survey from the SDSS (cf., Strauss 1997 for a description) will allow 
a determination of Pc{k), the galaxy power spectrum in redshift space. The quantity Pc{k) 
differs from the mass power spectrum P(k) due to two effects. First, the galaxy distribution 
may be biased with respect to the mass distribution. On large scales, models (Weinberg 
1995; Kauffmann, Nusser, & Steinmetz 1997; Scherrer & Weinberg 1997) indicate that the 
mass and galaxy density fields are directly proportional to one another. The proportionality 
constant is referred to as the galaxy bias parameter, b. Second, peculiar velocities cause 
the density contrast of galaxies in redshift space to appear appreciably stronger than in 
real space. In linear theory, the power spectrum is boosted by a factor 1 + |/3 + |/3 2 , 
where f3 = Q^/b (fl m is the matter density in units of the critical density of the Universe 
p c = 3Hq/(87tG), with H denoting the Hubble constant). The net result (Kaiser 1987; cf., 
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Hamilton 1998 for a review) is that on linear scales, the galaxy power spectrum is given by: 

Pa{k) = b 2 (l + 2 -(3 + i/? 2 ) P(k) = b 2 eff P(k). (3) 

The observables Cti and Pc{k) depend on the cosmological parameters Hq, (baryon 
density in units of p c ), Q m , (density contribution from the Cosmological Constant in 
units of p c ), and r ri (reionization optical depth) only through Atz(&, r = t ), b e jj, and T(k). 
The MAP data alone will allow rather tight constraints on the baryon/photon ratio (which 
determines Qbh 2 ), the matter /photon ratio (which determines Q m h 2 ), and the geometry of 
the universe, and test the basic inflationary scenario (Spergel 1994; Knox 1995; Jungman et 
al. 1996; Bond, Efstathiou & Tegmark 1997; Zaldarriaga, Spergel & Seljak 1997; Dodelson, 
Kinney & Kolb 1997; Lidsey et al. 1997; Copeland, Grivell & Liddle 1997; Efstathiou & 
Bond 1998; Eisenstein, Hu & Tegmark 1998). However, there is a near-degeneracy in several 
sets of parameters (Bond et al. 1994; Bond, Efstathiou Sz Tegmark 1997; Zaldarriaga, 
Spergel & Seljak 1997; |Huey fc Dave 19971) including the overall amplitude of the power 
spectrum at k ~ 0.1 /iMpc -1 and the optical depth. Since Ar;(fc, r = r ) and T{k) depend 
on the cosmological parameters differently, by combining the MAP and SDSS data, we can 
break these degeneracies and determine cosmological parameters to impressive accuracies. 
Several other recent papers have also considered how well cosmological parameters can be 
constrained with the combination of future CMB and galaxy survey data: Tegmark (1997) 
has studied a somewhat different and smaller set of the cosmological parameters than we 
consider here; Hu, Eisenstein, & Tegmark (1997) and Eisenstein, Hu, & Tegmark (1998) 
have tested the sensitivity of the data to the neutrino mass, which we do not consider in 
this paper. See also Lineweaver & Barbosa (1998) and Webster et al. (1998), who fit for 
cosmological parameters using the best available CMB and redshift data available. 

Both Cti and Pc{k) are proportional to the amplitude of the power spectrum of the 
primordial density fluctuations A 2 s (k). Because of the finite resolution of the MAP satellite, 
the errors on the Cti increase rapidly with wavenumber for I > 800, while the errors on 
Pc(k) from the SDSS decrease with wavenumber; MAP and SDSS thus quite naturally 
complement each other in the determination of A 2 s (k). Assuming that inflation occurred, 
all the obtainable information about the inflationary model is contained in A 2 s (k). 

Some inflationary models predict a significant gravity wave contribution to the 
microwave background fluctuation spectrum at small I. However, even in these models, 
gravity waves do not make a significant contribution to the microwave background 
fluctuation spectrum at I > 100 ( |Crittenden et al. 1993| , podelson, Knox, fc Kol5 



1991 Allen fc Koranda 1995| , |Wang 1996| , [Turner fc Wang 1996|) . Since most of the 



information on the cosmological parameters comes from relatively small angular scales 
(/ > 100), the inclusion of tensor perturbations would not change our results qualitatively. 

We parameterize A 2 s {k) by a series of steps equally spaced in log/c bins, each step with 
a constant amplitude in its bin. By taking these amplitudes to be independent parameters, 
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we probe features in the primordial power spectrum on all scales. The measurement of the 
primordial power spectrum independent of specific inflationary models from the combined 
MAP and SDSS data should shed light on our understanding of the early Universe. 

In Section 2, we present the statistical technique we use in estimating errors of 
parameters. Before considering a general primordial power spectrum, we determine the 
accuracy with which cosmological parameters will be measured with a pure power law 
A 2 s (k) in Section 3. In Section 4, we relax this constraint, and discuss the determination 
of the primordial power spectrum, using the stepwise parameterization described above. 
Section 5 contains a summary of the paper. 



2. The Fisher Information Matrix 

The Fisher information matrix of a given set of parameters, s, approximately quantifies 
the amount of information on s that we "expect" to get from our future data. The Fisher 
matrix can be written as 

F = -^± (4) 

where L is the likelihood function, the expected probability distribution of the observables 
given parameters s. The Cramer-Rao inequality ([Kendall fc Stuart 1969|) states that no 



unbiased method can measure the z-th parameter with standard deviation less than 1 / yRii 



if other parameters are known, and less than y(F _1 )jj if other parameters are estimated 
from the data as well. Note that the derivatives in Equation (§]) are calculated assuming 
that the cosmological parameters are given by an a priori model, and thus the errors on the 
parameters are somewhat dependent on the assumed model. 

The likelihood function L in principle can be that of the raw data; see, for example, the 
discussion in Tegmark et al. (1998) and Dodelson, Hui, & Jaffe (1997). But as Tegmark, 
Taylor, & Heavens (1997) describe, it is far more convenient to work with the data in a 
more compressed form. For the galaxy data, we work with the observed galaxy power 
spectrum, which can be determined using the techniques reviewed, e.g., in Strauss (1997) 
or Tegmark et al. (1998). The values of the measured power spectrum are approximately 
statistically independent if measured at values of k spaced at intervals of 2ir/R, where R 
is the characteristic length of the survey volume. Similarly, for a full-sky CMB mapping 
experiment like MAP, the individual C/'s are very nearly statistically independent (Oh, 
Spergel & Hinshaw 1998). 
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2.1. CMB Data 

MAP will measure the CMB polarization as well as temperature anisotropy. The 
CMB radiation field is described by a 2x2 intensity tensor 7^. The Stokes parameters are 
Q = (In — 722)/4 and U = Iu/2, while the temperature anisotropy is 5T/T = (In + -?22)/4. 
The temperature anisotropy can be expanded in standard spherical harmonics, while 
(Q ± iU) have to be expanded using spin weight ±2 spherical harmonics ±2^/m ( |Gelfand~et 



|al. 1963| , poldberg et al. 1967] , [Zaldarriaga fc Seljak 1997| ; see also |Kamionkowski, Kosowsky. 
fe Stebbins 1997 for an alternative expansion): 



8T/T(h) = ^a T ,im Y irn(n), 

I'm 

(Q + iU)(h) = ^a 2M2 Yi m (ri) 



1 1 n 

(Q-iU)(h) = J2a_2,im-2 Y irn(n), (5) 

I'm 

where fi is the unit direction vector in the sky. The even and odd parity linear combinations 

of a 2) i m and a_ 2 ,im are 

0'E,lm = — + 0-2,Zm)/2, 

0>B,lm = i(a2,lm ~ 0-2,Zm)/2. (6) 

Note that a E j m and a,B,im behave analogously to the electric and magnetic fields respectively 
under parity transformation. The statistics of the CMB are specified by the power spectra 
of the variables (a-zym, a E j m , a>B,im) together with their cross-correlations. For density 
(scalar) perturbations, the odd parity mode B vanishes. The power spectra can be written 
as flSeljak fc Zaldarriaga 19961 ) 

Cti = (\a TM \ 2 ) = (4n) 2 J^A 2 (k)\A Tl (k,r = r )\ 2 , 

C El = (\a EM \ 2 ) = (An) 2 J^Al(k)\A El (k,T = T )\ 2 , 

Cci = (a TM a EM ) = (4tt) 2 J — A 2 s (k) A Tl (k, r = r ) A El (k, r = r ), (7) 

where Ax, E i(k, r = tq) are integrals over conformal time r of the sources which generate the 
CMB temperature fluctuations and polarization. 

The assumption of random-phase density fluctuations does not necessarily imply that 
the errors on the Ct, e ,ci and the P(k q ) are Gaussian-distributed, although the Central 
Limit Theorem allows us to make this approximation. Following Jungman et al. (1996), 
Goldberg & Strauss (1998), Tegmark (1997), and Zaldarriaga, Spergel, & Seljak (1997), we 
can then write the Fisher matrix for the combined MAP and SDSS data as: 

p _ dCxi Coy- 1 (Cxi C±i) dCYl hV 1 dP G{k q ) dP G (k q ) ^ 
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where X, Y — T,E,C. Cov^ 1 (Cxi, Cyi) is the inverse of the covariance matrix 
Coy(Cxi, Cyi) = (ACxiACyi) ■ For each I, one has to invert the covariance matrix and 
sum over X and Y. <jp q are the standard errors in the measurement of P(k q ) (equation |12| 
below). We compute the Cxi and P(k) using the CMBFAST Boltzmann code by Uros 
Seljak and Matias Zaldarriaga flSeljak fc Zaldarriaga 1996] ). We calculate those derivatives 
which cannot be found analytically using two-sided finite differences. We take step sizes of 
2% of its model value for h, 5% of the model value of Q m for Q^, and 5% of the model value 
of each parameter for the other parameters. 

The covariance matrix of CMB power spectrum estimators has diagonal elements given 
approximately by ( [Seljak 199C , Zaldarriaga fc Seljak 1997 , Kamionkowski, Kosowsky, fc 
IStebbins 19971) : 



Cov(C^) 
Cov(C 2 El ) 



(C Ti + w T 1 B l 



Cov(C, 



ci) 



(21 + I) f sk y 

1 



c 2 



(21 + I) f sky 

while the off-diagonal matrix elements are: 



C 



TI 



w T 1 B l 



C 



El 



Wp X Bl 



(9) 



Cov(C T iC E i) 
Cov(CtiC C i) 
Cov(C El C C i) 



(21 + l)f sky 
2 

(21 + l)fsky 

2 



C 2 



(21 + 1)1 



sky 



Cci (C T i + w T 1 B l 
Cci (C E i + Wp^Bl 



(10) 



where f s k y is the fraction of sky which is mapped. We take f s k y = 0.8 for MAP. We have 
defined w7 T ^ = o\t p ffiwhm^ where ot.p is the noise per pixel in the temperature and the 
polarization measurements, and 6f w hm is the full width at half maximum of the beam in 
radians. The quantity Bf = exp — (0.425 6f whm I) 2 is the window function of the Gaussian 
beam. 



For MAP, both temperature and polarization data are obtained from the same 
experiment by adding and differencing the two polarization states, hence a\ = crp/2, and 
the noise and polarization signals are uncorrelated. MAP's resolution at its three highest 
frequency channels (40, 60, and 90 GHz) is anticipated to be 9f w h m = 0.47°, 0.35°, and 0.21° 
respectively, with the corresponding <tt=27 /zK, 35 /iK, and 35 fiK for its two year nominal 
mission (Bennett et al. 1997; prttp : / /map . gsf c . nasa . gov| ) . We do a weighted average in 
the three channels c, yielding 



Wt 



.2 Q2 

T,c fwhm,c 



(0.10 fiK) 



2 ' 
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w T /2, 

E c w c T Bl E c e- l2 <c/(a 2 T J 2 fwhmtC ) 

^ c t Ecl/(4A ta ,c) 

In our analysis, we ignore foreground contamination. Current estimates (Efstathiou & 
Bond 1998; Refregier et al. 1998) suggest that foregrounds will not significantly affect MAP 
temperature data. Foreground contamination, however, may significantly limit analyses of 
polarization data, particularly at low I (Keating et al. 1998). Thus, we report our results 
both for analyses that use only the temperature data and for analyses that use both the 
temperature and polarization data. 



Wp = 

Bf = 



2.2. Galaxy Data 

We take the SDSS standard error in the measured galaxy power spectrum from 
Feldman, Kaiser, & Peacock (1994). For the power spectrum determined over a volume in 
fc-space of 14: 

a P _ f W 3 /^^(r)^r)[l + =^] 2 y /2 



P G (k) I 14[/rf 3 rn 2 (r)^ 2 (r)] 2 j 



(12) 



where n(r) is the selection function, and ip{ r ) is the weighting function. For spherical 
shells, Vk q = ^(kq — fcg_i), where again, the spacing between the k g is 2n/R. We choose 
R = 1466/i- x Mpc for the SDSS, following Goldberg & Strauss (1998). We use the SDSS 
selection function calculated by Goldberg & Strauss (1998) using a mock SDSS redshift 
sample created by D. Weinberg (cf., Cole et al. 1998). We take the weighting function 
ip(r) = 1 throughout. Equation flT2| ) assumes that the power estimates in each k band are 
uncorrelated. Goldberg & Strauss (1998) show that this assumption causes at most a 10% 
underestimation of the true errors of parameters. 



3. Determination of the Cosmological Parameters 

In this section, we discuss how well the cosmological parameters can be determined 
by combining the MAP and SDSS data. Although we consider a general primordial power 
spectrum in Sec. 4, we here restrict ourselves to the simpler case of a power-law model 
for the primordial power spectrum generated by inflation, k A 2 s (k) oc k ns . We consider h, 

Qb, r r i, ris, Q m , C*2 = |«T,2m| 2 (which sets the CMB normalization), and b e g as free 
parameters. In Sec. 4, we will see how our determination of these quantities worsens once 
we give up a priori constraints on A 2 s (k). 



The Fisher matrix depends on the true underlying power spectrum. We consider here 
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three models, standard cold dark matter model (SCDM), the CDM model with cosmological 
constant (ACDM), and open CDM (OCDM). The model values of the cosmological 
parameters are shown in Table 1. The Ci values in Table 1 are obtained by normalizing the 
models to the COBE four- year data, following Punn fc White (1997)| . We take the galaxy 



bias parameter to be b = l/o"s, with as derived using the model C*2 value, and using the 
standard value for the observed fluctuations of galaxies within 8 /i _1 Mpc spheres ( Pavis & 



Peebles 1983[). Note that as we vary parameters to compute the Fisher matrix, we keep 



Vt m + Q\ = 1 for the SCDM and ACDM models, because these models are flat by definition. 
In the open CDM model, we allow for the simultaneous determination of Vt m and Q\. 

Table 1: Model values of cosmological parameters 





h 


Q A 


V b 


7~ri 


n s 


Vl m 


c 2 




b = l/a 8 


SCDM 


0.5 





0.05 


0.05 


1 


1 


1.139 x 10" 


-10 


0.83 


ACDM 


0.65 


0.7 


0.06 


0.1 


1 


0.3 


1.455 x 10" 


-10 


1.29 


OCDM 


0.65 





0.06 


0.05 


1 


0.4 


1.607 x 10" 


-10 


1.72 



In the linear theory of mass density fluctuations, the bias factor b is degenerate with 
the CMB normalization G^. Nonlinear effects in the present-day galaxy power spectrum 
break this degeneracy. We include nonlinear effects following Peacock & Dodds (1996). 
Given a linear power spectrum Pl(/cl), the effective local power law index can be defined as 



n e ff{k L ) = n L {ki 



d\nP L 



dliak 



k=k L /2 



The nonlinear power spectrum P/vl(/c) is related to the linear power spectrum by 

A 2 NL (k NL ) = f NL (A 2 L (k L ),tt), 
where A|(fc L ) = klP L (k L )/(2n 2 ), and = k 3 P NL (k)/(2n 2 ). 

1 + B-fx + [Ax] a ^ 



f NL (x,£l) = X 



([Ax} a g*(n)/[VxV 2 ]y 



l/ 7 



with 



(13) 



(14) 



(15) 



A = 


0.482 (l + n eff /3Y 


-0.947 


B = 


0.226 (1 + n eff /3Y 


-1.778 


a = 


3.310 (1 + n eff /3)" 


-0.244 


7 = 


0.862 (1 + n eff /3Y 


-0.287 


V = 


11.55(1 + 71^/3)" 


-0.423 



(16) 
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and ( Carroll, Press, fc Turner 1992 ) 



(7(fi) = ^ m [^ 7 -fi A +(l + fi m /2)(l + fi A /70)] \ (17) 

The nonlinear wavenumber is related to the linear wavenumber by 

k NL = k L [l + A 2 NL {k NL )} 1 '\ (18) 

Note that although our consideration of nonlinear effects is adequate in the context of this 
paper, it is inaccurate in several ways. First, even though we include nonlinear effects 
in P(k), we still use Kaiser's linear relation for redshift effects. Proper consideration of 
redshift distortions on small scales ( |Fisher fc JNusser 1996] , [Hatton fc Cole 1998| ) will lead 



to an effective bias b e g which depends on k as well as the cosmological parameters b and 
fl m . Second, we have explicitly assumed that the true bias b in real space is exactly linear 
and constant with scale; there must be scale-dependence and non-linearities, especially at 
small scales, where the density fluctuations are large (cf., Dekel & Lahav 1998; Blanton et 
al. 1998). Third, as the strength of the baryon-induced wiggles in the power spectrum is 
proportional to Qb/Q m , they are more pronounced in the ACDM and OCDM models. The 
Peacock & Dodds (1996) non-linear formulae only apply to smooth linear power spectra; 
we therefore smooth the linear power spectra of the ACDM and OCDM models (using 
the analytical formulae in |Eisenstein fc Hu 1998j) before applying the non-linear formulae. 



Since we are discarding the information contained in the baryon-induced wiggles, we obtain 
overestimates of the true errors. Finally, non-linear effects will cause mode-mixing between 
estimates of the power spectrum we have assumed here to be independent (cf., Jain & 
Bertschinger 1994). These effects are difficult to include in the Fisher matrix formalism, 
and we do not attempt to do so here. We carry out all calculations in this section to two 
limits in the k q of the galaxy power spectrum, k max = 0.1 /iMpc -1 and k max = 0.5 /iMpc -1 . 
The Peacock & Dodds (1996) formalism indicates that for the models considered here, 
non-linear effects become significant for values of k somewhat larger than 0.1 /iMpc -1 . So 
the former value of k max corresponds to a scale where non-linear effects are negligible, while 
the latter corresponds to a scale 2 7c/k « 12/i~ 1 Mpc where nonlinear effects are becoming 
strong. It is work for the future to check the validity of our calculations in this non-linear 
regime. 

To calculate the Fisher matrix, we need the derivatives of the Cxz's (X denotes T, E, 
C) and Pcik) = b^PMhik) with respect to the cosmological parameters. The derivatives 
with respect to Ci and b e g can be found analytically: 

dC xl 



dlnC 2 
dC xl 



C 



XI, 



d In b 



0. 



dP G (k) d\nf NL (x,Q) 
din Co " G[ ' 8\nx 
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dP G (k) 



2P G (k), 



d\nb eff 

where x = A 2 L {ki), and k^ is related to k via equation 
both Qa and Q m to be free parameters. Hence 



(19) 



In the OCDM model, we take 



dP G (k) 

on a 



P G {k) 



d In f NL (x,tt) 



+ 



d In f NL (x, ft) 



d In x 



Ma 



d In P L (k L ) 



(20) 



where Pi^l) depends on Q\ only through an overall normalization for the k range relevant 
to redshift surveys (on superhorizon scales, Pl(/cl) has a fc-dependent ^a dependence in 
the OCDM model). The derivatives with respect to the other cosmological parameters are 
calculated numerically by finite differences. 

Tables 2, 3, and 4 show the Fisher Matrix \-a error bars for the cosmological 
parameters, Asj = 



(F-% for SCDM, ACDM, and OCDM respectively. For each 
parameter, we list the errors for six cases: the MAP temperature data taken alone, the 
MAP temperature data with SDSS data to k max = 0.1/iMpc -1 , the MAP temperature 
data with SDSS data to k max = 0.5/iMpc -1 , and the corresponding results when MAP 
polarization data are included as well. 



Table 2: Fisher Matrix l-a error bars for the cosmological parameters of SCDM. 



1-0" error 


MAP T 


MAP T + 
SDSS 1 


MAP T + 
SDSS 2 


MAP( T+p ) 


MAP^+ 
SDSS 1 


MAP^ + 
SDSS 2 


AlnC 2 


0.17 


0.097 


0.016 


0.065 


0.053 


0.016 


A In ft 


0.052 


0.041 


0.021 


0.033 


0.028 


0.018 


AOv 


0.15 


0.11 


0.057 


0.091 


0.076 


0.047 


A\n(n b h 2 ) 


0.028 


0.024 


0.019 


0.020 


0.019 


0.017 


A In r ri 


2.4 


1.5 


0.46 


0.39 


0.38 


0.30 


A Inns 


0.017 


0.014 


0.0081 


0.0085 


0.0076 


0.0059 


Alub e ff 




0.064 


0.028 




0.063 


0.022 



Note. - - "T" denotes temperature only, "(T+P)" denotes temperature plus polarization 
for MAP data. "1" and "2" denote the SDSS cutoff k max = 0.1/iMpc" 1 , and 0.5/iMpc- 1 
respectively. 



The peak features (called "Doppler peaks" , "acoustic peaks" , or "Sakharov peaks" ) in 
the CMB angular power spectrum Cti are due to acoustic oscillations in the baryon-photon 
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Table 3: Fisher Matrix l-o error bars for the cosmological parameters of ACDM. 



1-cr error 


MAP T 


MAP r + 
SDSS 1 


MAP r + 
SDSS 2 


MAP( T+p ) 


MAP™+ 
SDSS 1 


MAP™+ 
SDSS 2 


AlnC 2 


0.21 


0.15 


0.010 


0.043 


0.043 


0.010 


Aln/i 


0.066 


0.048 


0.025 


0.032 


0.030 


0.021 


Alnfi A 


0.076 


0.055 


0.027 


0.037 


0.035 


0.023 


A\n(n b h 2 ) 


0.044 


0.034 


0.022 


0.021 


0.021 


0.018 


A 111 T r i 


1.33 


0.95 


0.18 


0.18 


0.18 


0.13 


A In ns 


0.035 


0.026 


0.014 


0.014 


0.013 


0.011 


A\nb eff 




0.076 


0.022 




0.038 


0.021 



Table 4: Fisher Matrix 1-cr error bars for the cosmological parameters of OCDM. 



l-a error 


MAP T 


MAP r + 
SDSS 1 


MAP T + 
SDSS 2 


MAP( T+p ) 


MAP™+ 
SDSS 1 


MAP™+ 
SDSS 2 


AlnC 2 


0.13 


0.10 


0.034 


0.081 


0.064 


0.033 


Aln/i 


0.11 


0.085 


0.018 


0.093 


0.069 


0.015 


A \nfl m 


0.13 


0.099 


0.011 


0.13 


0.083 


0.010 


AOv 


0.088 


0.060 


0.010 


0.085 


0.055 


0.0097 


A\n(tt b h 2 ) 


0.041 


0.033 


0.017 


0.038 


0.031 


0.017 


A In T ri 


1.5 


1.5 


0.46 


0.57 


0.57 


0.37 


A Inns 


0.026 


0.021 


0.0079 


0.023 


0.020 


0.0076 






0.14 


0.029 




0.14 


0.026 



fluid at the time of last scattering of the CMB photons. As mentioned above, the amplitude 
of these peaks is determined by VL^h 2 and VL m h 2 , while their location is determined by the 
geometry of the universe (roughly Q m + Q\ in a nearly flat universe), hence there is a near 
degeneracy in one dimension of the parameter space of h, fib, and fl m (Bond et al. 1994; 
Zaldarriaga, Spergel & Seljak 1997; Huey fc Dave 1997 ). For models with only six or seven 



free parameters, this degeneracy is not severe (see Table 2); however, it is significant if we 
consider models with many more free parameters (e.g., a variable spectral index, or massive 
neutrinos) . 

Large scale structure observations are sensitive to an orthogonal set of parameters. 
The turnover in the matter power spectrum is set by the horizon size at matter-radiation 
equality (e.g., Peebles 1993). Thus, SDSS will make an accurate measurement of Q m h. 
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The Silk damping scale, and the acoustic oscillations also affect the transfer function, so 
that VL b h can also be determined from the SDSS observations (Goldberg & Strauss 1998). 
Thus combining the MAP and SDSS data allows us to do two important things. First, to 
the extent that parameters derived from the two independent data sets are consistent, we 
obtain a powerful check of the correctness of the physical model which ties the two together: 
gravitational instability theory, the assumption that the fluctuations are adiabatic, and the 
linear biasing paradigm. This is illustrated in Figure 1, which shows the 68.3% and 95.4% 
confidence contours in the Q^-h (a) and Q m -h (b) planes. The dotted lines are for MAP 
temperature data only, the dashed lines are for SDSS data only [k max = 0.5/iMpc -1 ), and 
the solid lines are for the combined MAP temperature and SDSS data. The SCDM model 
is assumed. Second, the joint analysis of the two independent data sets breaks some of 
the degeneracies inherent in the MAP data taken alone and improves the accuracy in the 
determination of parameters. Figure 2 shows the 68.3% and 95.4% confidence contours in 
the Q\-Q m plane in the OCDM model, with the same line types as in Figure 1. Note that 
using MAP data alone, there is degeneracy between Q m and Q\; using SDSS data alone, Q\ 
can not be determined because it is degenerate with C*2 (see equation EDI). The combined 



data give a much better determination of Q m and than does either the MAP data or the 
SDSS data alone. 

The CMB temperature spectrum is also approximately degenerate in its dependence 
on the reionization optical depth r r j and the CMB normalization C*2, because reionization 
suppresses the spectrum by a factor of exp(— 2r ri ) on small scales relative to the large 
scales. Figure 3 shows the confidence contours in the r ri -C2 plane, with the same line types 
as in Figure 1. Combining MAP data with SDSS data dramatically breaks the degeneracy 
between r r j and C*2 by adding information on the small scales. 

Figures 1-3 show MAP temperature data only, because the extent to which MAP 
polarization data can be used may be limited by foregrounds (see Section 5). Also, if we 
allow a cutoff of k max = 0.5/iMpc" 1 (as we do for the rest of this paper), adding MAP 
polarization data to the combined MAP and SDSS data does not have a significant impact 
on the accuracy of the determination of cosmological parameters (see Tables 2-4) . 

Note that the errors of the parameters from combined MAP and SDSS data are very 
sensitive to the SDSS cutoff k max . This is as expected, since the small scale information 
from the SDSS is determined by its cutoff k max . Indeed, in most cases, the galaxy data to 
kmax = 0.1 /iMpc -1 gives only a moderate improvement to the constraints the MAP data 
alone give; much of the improvement comes from the non-linear galaxy regime. 

To summarize, when MAP temperature data is considered alone, there is strong 
degeneracy between the overall amplitude of the matter power spectrum C<i and the 
reionization optical depth r r j, and between the matter density fraction Q m and the density 
fraction contribution from cosmological constant Q\. Combining MAP and SDSS data 
dramatically breaks the degeneracy between C 2 and r ri , and reduces the degeneracy 
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between Vt m and f^A, leading to a significant improvement on the determination of C2, r ri , 
Q m , and f^A, enabling us to measure the effective bias between the matter power spectrum 
and the galaxy redshift power spectrum b e ff. The accuracy in the determination of other 
cosmological parameters is also improved by combining MAP and SDSS data, but to a less 
impressive extent. In the next section, we study the measurement of the primordial power 
spectrum itself, which can only be achieved by using the combined MAP and SDSS data. 



4. Measurement of the Primordial Power Spectrum 

In the previous section, we considered a power-law primordial power spectrum, and 
derived the constraints on cosmological models that the combination of MAP and SDSS 
data will give. We now consider a more general primordial power spectrum; as we discussed 
in the Introduction, various inflationary models make a range of predictions for the form of 
the primordial power spectrum. We parameterize the primordial power spectrum A 2 s (k) as 

A 2 s(k) = l aU f X t <k \ t ■ , (21) 
&w [ Oj, for < k < fy, 1 > 1 

where k± = 0.001 /iMpc -1 , and ki (i = 2,3, ...,20) are equally spaced in log/c from k\ to 
k = 0.5/iMpc -1 . The a, (i = 1, 2, 20) are taken to be independent variables. To compute 
errors in the determination of the a«, we need to know their true values, as well as true 
values of the cosmological parameters. We first consider the case that the cosmological 
parameters are known a priori; we consider the realistic case in which we know neither 
the primordial power spectrum nor the cosmological parameters later in this section. For 
definiteness, we adopt a model in which the true primordial power spectrum is A 2 s (k) = 1, 
and the cosmological parameters are those of SCDM (see Table 1). 

In the Fisher matrix approach we have adopted in this paper, we need to compute the 
derivatives of the observables with respect to the parameters, which now include the values 
of the steps a, which make up the primordial power spectrum. The derivatives of the Q's 
with respect to cij are taken to be finite differences with step-size of 5% the model value of 
aj. The derivatives of the galaxy power spectrum with respect to eij are given by 



dlnP G (k) = f dln a f » L J x) ±, for k < h (i = l),or k^ <k<k i {i> 1), 
(9a,' I 0, otherwise, 



(22) 



where x and /atl(x) are the same as in Equation (|T^). 

In this section, we take the cutoff wavenumber of the SDSS to be k = 0.5Mpc _1 
throughout. Also, when we consider the combined MAP and SDSS data, only the 
temperature data from MAP is used; adding the polarization data does not have a 
significant effect on errors from the combined MAP and SDSS data, while it is of great 
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interest to consider the SDSS data as an alternative to the MAP polarization data in 
constraining inflationary models. 

Figure 4 shows the accuracy in the determination of the bin amplitudes of the primordial 
power spectrum, k A 2 s (k); the 1-a error bars are shown for MAP temperature data only (a), 
SDSS data only (b), and MAP temperature and SDSS data combined (c). Assuming the 
cosmological parameters are known, for 20 bins in log k (k < 0.5/iMpc -1 ), the primordial 
power spectrum can be determined to around 16% accuracy for k ~ 0.01 /iMpc" 1 , and to 
around 1% accuracy for k ~ 0.1 /iMpc -1 (see Figure 4(c)) by combining MAP temperature 
and SDSS data. 

Not surprisingly, the determination of the primordial power spectrum bin amplitudes 
does depend on the a priori knowledge of the cosmological parameters. We now consider the 
realistic case in which the cosmological parameters are not known a priori; in particular, in 
addition to the a iy we also solve for the four parameters h, Q\, and r ri . The primordial 
power spectrum bin amplitudes and the cosmological parameters can not be simultaneously 
determined using MAP temperature data only, because r ri is degenerate with the for 
k > 0.01 /iMpc -1 . However, polarization data removes this degeneracy. Figure 5(a) shows 
the resulting error bars on the a« if both the temperature and polarization data from MAP 
are used; they are significantly larger than in Figure 4(a) where only temperature data is 
used. 

Figure 5(b) shows the effect of combining the MAP temperature data with the SDSS 
data in the analysis, adding the galaxy bias b e g as a fifth parameter. The uncertainty in 
the determination of the primordial power spectrum bin amplitudes aj increases by a factor 
up to 3 on small scales relative to the case in which the cosmological parameters are known 
a priori (Figure 4(c)). This difference is indistinguishable in the figures, because it enters 
through error bars that are only of the order of a few percent. The errors in the are 
not significantly increased by the inclusion of cosmological parameters in the parameter 
estimation, because most of the information on the cosmological parameters comes from 
small angular scales, where the power spectrum can be determined to great accuracy 
from the SDSS data. Comparison of Figure 5(a) and Figure 5(b) illustrates the critical 
importance of combining the MAP data with the SDSS data in the determination of the 
primordial power spectrum bin amplitudes. 

Note that our choice of 20 bins in log/c to parameterize A 2 s (k) is somewhat arbitrary. 
Obviously, the errors on the bin amplitudes would increase if we were to increase the 
number of bins. When the cosmological parameters are estimated together with the A 2 s (k) 
bin amplitudes a iy we expect the covariance between the cosmological parameters and the 
Oj's to increase rapidly with smaller bin size, as the primordial power spectrum becomes 
flexible enough to mimic the Doppler peaks and other features of the C\. 

To illustrate how the errors on the cosmological parameters change as they are 
estimated jointly with primordial power spectrum bin amplitudes using the combined MAP 
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temperature and SDSS data, we show in Table 5 the Fisher matrix 1-er error bars for two 
groups of cosmological parameters, (h, f^, T ri, b e g) and (h, Q\, Oj,, r ri , b e jj), when they 
are estimated jointly within each group and when each group is simultaneously estimated 
with 20 independent bin amplitudes of the primordial power spectrum. The cosmological 

Table 5: Fisher Matrix l-a error bars for the cosmological parameters. 



l-a error 


4 parameters 


4 parameters & 20 bin 
amplitudes of A 2 s {k) 


5 parameters 


5 parameters & 20 bin 
amplitudes of A 2 s {k) 


Aln/i 


0.0048 


0.0062 


0.017 


0.025 








0.046 


0.064 


A\n(n b h 2 ) 


0.016 


0.039 


0.016 


0.046 


A In T r i 


0.056 


0.19 


0.12 


0.27 


A\nb eff 


0.0052 


0.0094 


0.012 


0.017 



parameters can still be determined to impressive accuracy even when they are estimated 
simultaneously with 20 independent bin amplitudes of the primordial power spectrum! 
Of course, we do not include the overall amplitude of the power spectrum C2 as a free 
parameter, as it is degenerate with the by definition. 

Figures 4 and 5 assume that the true primordial power spectrum is A 2 s (k) = 1. For the 
purpose of illustration, we can choose our model primordial power spectrum A 2 s (k) to have 
broken scale-invariance based on the Holman et al. (1991) and Adams et al. (1997) models. 
We write 

f 1, k < 0.01 /iMpc -1 , 
A 2 s (k) = l A + Bk } 0.01 /iMpc -1 < k < 0.1 /iMpc -1 , (23) 
{ 0.25 AT ' 2 , k > 0.1/iMpc" 1 , 

where A and B are chosen such that A 2 s (k) is continuous. Figure 6 shows the primordial 
power spectrum of Equation ([23|) with l-a error bars for the combined MAP temperature 
and SDSS data. Clearly, the accuracy in the determination of A 2 s {k) is not sensitive to the 
assumed underlying model; the errors on A 2 s (k) bin amplitudes do not change significantly 
with model where the errors are larger than of order a few percent. 

There have been recent claims of possible observational detection of features in the 
matter power spectrum on the scale of k ~ 0.1 /iMpc" 1 ( [Peacock 19371 ). Combining MAP 



and SDSS data will enable us to confirm such features, and to determine whether they are 
intrinsic to the primordial power spectrum, independent of specific inflationary models. 
Indeed, it would be straightforward to translate our bin amplitude errors of the primordial 
power spectrum into constraints on existing inflationary models; our errors are small 
enough to be very interesting. However, it is perhaps more exciting to anticipate that the 
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measurement of the primordial power spectrum in wavenumber bins from the combined 
MAP and SDSS data will provide a lead to the correct inflationary model. 



5. Summary 

Cosmic microwave background experiments measure fluctuations in the curvature 
of space at the surface of last scattering. If the primordial fluctuations are adiabatic, 
then MAP will be a sensitive probe of the geometry of the universe, the matter/photon 
ratio (Q m h 2 ), the baryon/photon ratio (Qbh 2 ) and the primordial power spectrum for 
k < 0.2/iMpc- 1 . 

Observations of large scale structure measure the large scale distribution of galaxies. 
If the primordial fluctuations are adiabatic, then SDSS will be a sensitive probe of the 
primordial power spectrum for k > 0.02 /iMpc -1 and Q m h, which determines the horizon 
size at matter-radiation equality and thus the shape of the processed power spectrum. 
SDSS will also be a powerful probe of the nature of the dark matter (Hu et al. 1997). 
Because of the large number of galaxies with detailed photometric and spectroscopic data 
in the SDSS, we can determine the power spectrum and bias factor for several independent 
samples of galaxies, defined by morphology, color, or spectral characteristics. To the extent 
to which these independent samples yield consistent power spectra, we can determine the 
range in k over which the linear bias model is valid. 

The combination of the SDSS and the MAP data will be a powerful test of our basic 
cosmological models. Both experiments will accurately determine the amplitude of density 
fluctuations at 0.2 > k > 0.02 /iMpc -1 . Will there be a set of cosmological parameters and 
primordial power spectrum that is consistent with both sets of observations? If so, then we 
will have tested the gravitational structure formation paradigm, our interpretation of the 
primordial fluctuations as adiabatic, and the linear biasing paradigm. If not, this conflict 
will likely lead us to a deeper understanding of the origin of structure in the universe. 

In particular, from the combination of the MAP and SDSS data, we can obtain 
a measurement of the primordial power spectrum without reference to any specific 
inflationary model, assuming that the primordial fluctuations are adiabatic. If the values 
of the cosmological parameters that best fit the data are close to what we expect, but 
the primordial power spectrum differs significantly from the predictions of the current 
inflationary models, then it will be an indication of new physics in the early universe, and 
it will provide a solid starting point for building new inflationary models. 

If there is a theory consistent with both data sets, then the combination of the two 
observations will break the degeneracies inherent in each individual observation. Cosmic 
microwave background observations by themselves measure a combination of the amplitude 
of the primordial power spectrum and the optical depth. Large scale structure observations 
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measure the product of the current matter power spectrum and the bias parameter. In 
linear theory, the offset between the power spectrum as determined by the two sets of 
observations determines a combination of the bias parameter and the optical depth of the 
universe. Because of non-linear effects on the spectrum, this offset is scale-dependent, this 
will enable us to independently measure the bias parameter and the optical depth of the 
universe. The combination of the two experiments will also improve our ability to determine 
other cosmological parameters, in particular, the matter density and the cosmological 
constant. 

Measurements of the cosmic microwave background polarization can provide an 
independent measure of the optical depth of the universe r r j (Zaldarriaga, Spergel & Seljak 
1997). However, emission from both polarized galactic dust and synchrotron emission 
(Keating et al. 1998) may swamp the primordial polarization fluctuations at large angular 
scales. These foregrounds may significantly limit our ability to extract useful cosmological 
information from polarization measurements. Hence it is important that we can use the 
SDSS data as an alternative and independent aid to MAP in the determination of r ri . 

The Sloan Digital Sky Survey, and other large scale structure surveys, will make other 
independent measurements that will probe cosmological parameters and the primordial 
spectrum: redshift distortions, cluster properties, small scale velocity fields, and the 
evolution of structure. These will provide additional tests of the basic model and will 
further enhance our ability to measure cosmological parameters. 

With this in mind, there are a number of improvements that could be done in this 
analysis. The redshift distortions could be measured directly, and the effects could be 
included directly into the Fisher matrix analysis (cf., de Laix & Starkman 1998 and Hatton 
& Cole 1998 for a discussion of how well redshift-space distortions can be measured from 
the SDSS data). Similarly, we could parameterize the effects of galaxy and clustering 
evolution, and include these as parameters in the analysis. More challenging will be a 
proper accounting of non-linear effects on small scales, in the power spectrum, the redshift 
space distortions, and the bias (here assumed linear and independent of scale). The 
analyses in Tables 2-4 show that going to k max = 0.5 /iMpc" 1 gives us particularly strong 
constraints on cosmological parameters, but we have taken the non-linear effects into 
account in only a relatively crude way in this paper. It remains to be seen the extent to 
which the uncertainties due to these inaccuracies strongly affect the error bars we derive on 
the primordial power spectrum. 

The next few years will be a very exciting time in cosmology. 
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Fig. 1. — The confidence contours (68.3% and 95.4%) in the Q b -h (a) and Q m -h (b) planes 
assuming the SCDM model. The dotted lines are for MAP temperature data only, the dashed 
lines are for SDSS data to k max = 0.5 /iMpc -1 only, and the solid lines are for the combined 
MAP temperature and SDSS data. 

Fig. 2. — The confidence contours (68.3% and 95.4%) in the Q\-Q m plane for the OCDM 
model. The line types are the same as in Figure 1. 

Fig. 3. — The confidence contours (68.3% and 95.4%) in the T r i-Ci plane for the SCDM 
model. The line types are the same as in Figure 1. 

Fig. 4. — The primordial power spectrum with 1-cr error bars for MAP temperature data 
only (a), SDSS data only (b), and MAP temperature and SDSS data combined (c), assuming 
that the cosmological parameters are known to be that of the SCDM mo del (with A 2 s (k) = 1) 
listed in Table 1. The error bars in (b) which extend to the edges of the figure are infinite, as 
these refer to measurements of the power spectrum on scales larger than those which SDSS 
probes. 

Fig. 5. — The primordial power spectrum with 1-cr error bars, (a) Combined MAP 
temperature and polarization data, with four cosmological parameters (h, Q\, Qb, T r i) 
included in the parameter estimation; (b) Combined MAP temperature and SDSS data, 
with five cosmological parameters (h, fl^, th, b) included in the parameter estimation. 

Fig. 6. — The primordial power spectrum of Equation fl23f) with 1-cr error bars for the 
combined MAP temperature and SDSS data. 
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Fig. 1. — (a) The confidence contours (68.3% and 95.4%) in the VL^-h plane assuming the 
SCDM model. The dotted lines are for MAP temperature data only, the dashed lines are 
for SDSS data to k max = 0.5/iMpc _1 only, and the solid lines are for the combined MAP 
temperature and SDSS data. 
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Fig. 1. — (b) The confidence contours (68.3% and 95.4%) in the fl m -h plane. The line types 
are the same as in (a). 
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Fig. 2. — The confidence contours (68.3% and 95.4%) in the fl\-fl m plane for the OCDM 
model. The line types are the same as in Figure 1. 
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Fig. 3. — The confidence contours (68.3% and 95.4%) in the r ri -C 2 plane for the SCDM 
model. The line types are the same as in Figure 1. 



-30- 




Fig. 4.— The primordial power spectrum with 1-a error bars, assuming that the 
cosmological parameters are known to be that of the SCDM model (with A 2 s (k) = 1) listed 
in Table I. (a) MAP temperature data only. 
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Fig. 4. — (b) SDSS data only. The error bars which extend to the edges of the figure are 
infinite, as these refer to measurements of the power spectrum on scales larger than those 
which SDSS probes. 
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Fig. 4. — (c) MAP temperature and SDSS data combined. 
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Fig. 5.— The primordial power spectrum with 1-a error bars, (a) Combined MAP 
temperature and polarization data, with four cosmological parameters (h, Q^, ^b, T ri) 
included in the parameter estimation. 
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Fig. 5.— (b) Combined MAP temperature and SDSS data, with five cosmological 
parameters (h, Q\, Q&, r ri , b) included in the parameter estimation. 
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Fig. 6.— The primordial power spectrum of Equation ( p3|) with 1-a error bars for the 
combined MAP temperature and SDSS data. 



